Primary Data source: church_locations_panel.dta 

keep if church2=="Wesleyan Methodist Church" & year==1932
collapse (sum) exist, by(gridcell)
drop if exist==0 

keep if church2=="Presbyterian" & year==1932
collapse (sum) exist, by(gridcell)
drop if exist==0 

keep if church2=="Catholic" & year==1932
collapse (sum) exist, by(gridcell)
drop if exist==0 

The maps were created in QGIS by joining the dbf files to the <ghana_grid> shapefile.